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\0 • ABSTRACT 

I We consider the excitation of radial and non-radial oscillations in low-mass B stars by 

the iron-bump opacity mechanism. The results are significant for the interpretation 
of pulsations in subdwarf B stars, helium-rich subdwarfs and extreme helium stars, 
including the EC14026 and PG1716 variables. 

We demonstrate that, for radial oscillations, the driving mechanism becomes ef- 
fective by increasing the contrast between the iron-bump opacity and the opacity from 
. other sources. This can be achieved either by increasing the iron abundance or by de- 

' creasing the hydrogen abundance. The location of the iron-bump instability boundary 

I is found to depend on the mean molecular weight in the envelope and also on the ra- 

■ dial order of the oscillation. A bluer instability boundary is provided by increasing the 

iron abundance alone, rather than the entire metal component, and is required to ex- 
plain the observed EC14026 variables. A bluer instability boundary is also provided by 
higher radial order oscillations. Using data for observed and theoretical period ranges, 

O [ we show that the coolest EC14026 variables may vary in the fundamental radial mode, 

but the hottest variables must vary in modes of higher radial order. 
^ ' In considering non-radial oscillations, we demonstrate that g-modes of high radial 

' order and low spherical degree (l < 4) may be excited in some blue horizontal branch 

stars with near-normal composition {Z = 0.02). Additional iron enhancement extends 
the g-mode instability zone to higher effective temperatures and also creates a p-mode 
rS I instability zone. The latter is essentially contiguous with the radial instability zone. 

■ With sufficient iron, the p-mode and g-mode instability zones overlap, allowing a small 

region where the EC14026 and PG1716-type variability can be excited simultaneously. 
The overlap zone is principally a function of effective temperature and only weakly a 
function of luminosity However its location is roughly 5000K too low compared with 
the observed boundary between EC14026 and PG1716 variables. The discrepancy 
cannot be resolved by simply increasing the iron abundance. 

Key virords: stars: oscillations, stars: subdwarfs, stars: horizontal branch, stars: chem- 
ically peculiar, stars: early-type, stars: variables: other 
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1 INTRODUCTION 

Tlie discovery of a significant contribution to stellar opacity 
from iron-group element s at temperatures around 200 OOP K 
jRoeers fc Idesiaj Il992l: Ixiie Opacity Project Teaiiil Il995l) 
had major consequences for studies of pulsation in stars. 
First, it solved the long-standing Cepheid mass problem 
jMoskalik. Buchler fc Maromll99alKanbur fc Simonll994l) . 
Second, it provided a natural explanation for hitherto un- 



explained stella r variabilit y , notab ly with respect to the /3 
Cepheids (iDziemb owski fc Pamiatnvkh 1993) and some ex- 
treme helium stars Jsaicll^9^ . Third, it initiated searches 
for variability (predicted and observed) in stars hitherto 
thought to be stable. There are now known to be several 
classes of variable in which pulsations are driven by the 
opacity (or n—) mechanism excited by the Fe- or Z- opacity 
bump. It is one goal of this paper to explore some of the 
connections between these groups and hence to understand 
better the nature of Fe-bump driven pulsations. 



* E-mail: csj@arm.ac.uk 

f E-mail: saio@astr.tolioku.ac.jp 



An early success was the explanation of radial p ulsa- 
tions in the extreme helium star V652Her (ISaidll993) . to- 
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Figure 1. The loci of low-mass variable stars believed to pulsate 
due to Fe-bump instability, including EC 14026 variables (upright 
triangles), PG1716 variables (inverted triangles) and helium-rich 
variables (squares). A selection are labeled. HS0702-I-6843 and 
Balloon 090100001 exhibit both EC14026 and PG1716-type be- 
haviour iSchuh ct al. 2006; Orciro ct al. 2005). Straight diagonal 
lines represent logg = 6.0 (below) and 4.0 (above). These features 
are reproduced in Figs.|5|to|S| 



gether with the prediction llSai ol ll994l) and subsequent dis- 
covery of radial pulsatio ns in a second helium star, BXCir 
jKilkennv fc KoenI llQQsl) . In these stars, both with peri- 
ods of ~ 9 300s and with masses 0.59 and 0.47 }Ap) re- 
spect ively Jjeffery. Woolf fc Pollacccll200ll: IWoolf fc Jeffervl 
l2002r (Fig. 0, the cause of pulsation instability is clearly 
the Fe- opacity bump, exaggerated by a reduction of the 
background opacity due to hydrogen at these tempera- 
tures. The failure to detect any variability in a third helium 
star, HD144941, of similar temperature and luminosity but 
very low metallicity ij.Tefferv fc Hilll99fil : l,Tefferv fc HarrisonI 
Il997l) . confirmed the role of the Fe-opacity in driving pulsa- 
tions in these low-mass early-type stars. 

Another notable success has been the prediction and 
discov ery of short-period puls a tions i n subluminous B (sdB) 
stars llC harpinet et all 1 19961 Il997l : iKilkenriv et al.l Il997t 
iBilleres e t al. 199?}). Briefly, multi-periodic oscillations with 
periods between 90 and 600 seconds are observed in approx- 
imately 10% of hot subdwarf B stars; such stars are vari- 
ously known as EC14026 variables (after the class prototype 
EC 14026-2647) or sdBVs. These pulsations are success- 
fully explained if the stars are "extreme horizontal branch 
stars" , that is they consist of a helium-burning core of some 
0.5 M© overlaid by a very thin layer of hydrogen. Because 
of the high effective temperature and the high surface grav- 
ity, theoretical models have shown that radiative forces on 
the ions in the stellar envelope act differentially such that 
substantial chemical gradi ents are established over a dif- 
fusion time scale ~ lO^y ijMichaud et alJll98fl^ . The con- 
sequent levitation and accu mulati on of iron in layer s at 
around 200 000 K dChaver. Fontaine fc WesemaeJll993) en- 



hances the Fe-opacity bump sufficiently that radial and non- 
radial p-mode oscillations are excited. The theory success- 
fully explained the observed distribution (Fig. ^ of pul- 
sating and non-puls ating sdB stars in effect ive temperature 
and surface gravity dCharpinet et al]l200J) . Although con- 
sequential more on frequency eigenvalues than stability cri- 
teria, the theory is also sufficiently well developed that it 
has been possible to compare predicted and observed pulsa- 
tion frequencies in some puls ating sdB stars (Brassard et alJ 
bOOltlCharpinet et aljEoOSallbl) . 

The unexpected discovery of oscillations with periods of 
a few hours in many sdB stars lying red- ward of the EC14026 
instability domain presented a new challenge l|Green et alJ 
l2003lf) . While radiative levitation of iron is still operative in 
these stars, sometimes known as PG1716 variables (after the 
prototype PG 1716-1-426), p- modes are reported to be stable 
in the chemically stratified models (^haroinct ct al. 200lj). 
On the other hand, non-radial g-modes of high radial order 
(fe > 10) and high spherical degree {I > 3) were found to be 
unstable, but only in models of star s cooler than those i n 
which variability had been detected (|Fontaine et aljf2003l ). 
While the observed periods imply a g-mode origin, modes of 
such high degree are not generally thought to be observable 
as variations in total flux due to geometric cancellation. The 
challenge is therefore to shift the g-mode blue edge to higher 
effective temperatures and to excite modes of lower spherical 
degree. 

A further challenge has been provided by the detection 
of variability in the helium-rich sdB star L S IV— 14°1 16 with 
periods of ~ 1800 - 3000s (Ahmad & Jeffervl l200,'j ). While 
the canonical sdB stars mostly have surface helium abun- 
dances one tenth of the solar value (~ 10% by number), 
the much rarer He-rich sdB stars comprise a spectroscop- 
ically distinct and highly inhomogeneous group with sur- 
face helium abundances ranging from some 30% to nearly 
100% by number, and surface gravities ov er a much broader 
range than seen in normal sdB stars jAhmad fc Jeffervl 
[2003). The natural response to the discovery of pulsations 
in LSIV— 14°116 was to try to explain them in terms of 
Fe-bump instability in either a helium star like V652Her, 
or in some sort of mutant sdB star. With Teff similar to the 
EC14026 stars, but with a lower surface gravity (higher L/M 
ratio), the pulsation periods were too long to be exp lained 
by Fe-bump driven p- modes jAhmad fc Jeffervl l2005'). One 
question is whether g-modes might be excited in such a star. 
This challenge is important because a viable picture of the 
star that explains both the L/M ratio and the high surface 
abundance of helium has yet to be established. Any model 
that successfully reproduced the observed oscillations would 
assist this process. 

In order to address these questions, we need to develop 
our understanding of instability in low-mass stars. By gain- 
ing a general insight into what affects the type of pulsation 
that can be driven, it should be possible to determine what 
model properties need to be modified in order to reproduce 
observed oscillations. 

The important theme is clearly Fe-bump instability. 
We have already investigated the role of envelope hydro- 
gen abundance in t he excitation of Fe-bump pulsations 
jjefferv fc Saiol^l998^ . We recall that pulsational instability 
exists in all stars of sufficiently high L /M ratio due to the 
presence of strange-mode instability (iGautschv fc Glatzell 
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I1990I) . Fe-bump instability sets in for lower L/M stars with 
Tes ~ 20 — 30 OOOK when either the iron-abundance is 
raised sufficiently or the hydrogen abundance is reduced suf- 
ficiently. The second occurs because reducing the hydrogen 
opacity increases the contrast between the iron opacity and 
the background opacity due to other sources, effectively in- 
creasing the opacity gradients kt = {d\n K/dlnT)p, Kp = 
(9 In k/9 In p)t in the driving zone. 

Therefore there should be a physical connection be- 
tween, for example, the radial pulsations seen in the extreme 
helium star V652Her (essentially normal iron, reduced hy- 
drogen) and the p-mode oscillations in EC f 4026 variables 
(normal hydrogen, enhanced iron). Note that radial modes 
are special cases of p-modes with spherical degree I — 0; 
such modes are generally present in EC14026 stars alongside 
modes of higher spherical degree. However, the EC14026 in- 
stability region reported by Cha rpinct ct al. (2001.) lies con- 
siderably blue-ward of the Fe-b ump instability finger de- 
scribed bv lJefferv fc Saiol (Il998ll . Understanding this para- 
dox represents an important step towards understanding Fe- 
bump instability in general. 

Our initial goal is therefore to explore Fe-bump insta- 
bility for radial modes in comparatively simple (i.e. homo- 
geneous) stellar envelopes as a function of composition (sec- 
tion II). With success here, it becomes possible to restrict 
the volume of model space required to explore the stability 
of modes of higher degree, in particular g-modes (section 
III), and also to identify what is required in more detailed 
models in order to reproduce the observations accurately. 



2 RADIAL MODES 

Linear non-adiabatic models of radial pulsations have 
been constructed following the method described by ISaid 
il99d[). but including the most recent OPAL95 opaci- 
ties dlglesias &: Rogerjri99(j) . In all cases the equation of 
state (EOS) treats hydrogen, helium and a representative 
"metal", namely carbon. The contribution to the EOS of 
varying the heavy element distribution within "Z" is negli- 
gible at temperatures above 10 000 K since both the electron 
density and the mean atomic weight are completely domi- 
nated by hydrogen and helium. 

The parameter range considered is M/ Mq — 0.4 — 0.9, 
logTcft/A' = 4.20 - 4.60 and log L/M = 0.4 - 3.8 (solar 
units). In terms of composition, the hydrogen mass frac- 
tion X = 0.0 — 0.9 includes extremely helium-rich, normal 
and helium-depleted mixtures. The treatment of metal mass 
fraction Z is discussed below. The choice of L/M rather 
than L as independent variable was made because pulsation 
properties scale roughly with M (actually M^^^), simplify- 
ing the generation of model grids and the inter-comparison 
of results for different masses. 

For each combination of M,Tc{!,L/M,X and Z, a ho- 
mogeneous static stellar envelope is integrated assuming a 
standard surface boundary condition. This integration is ter- 
minated when the total pressure exceeds 10^^ dyn or the 
mass variable < 0.6M, whichever occurs first. This limit 
is justified for a stability analysis of radial modes (and p- 
modes in general) since amplitudes decay exponentially to- 
wards the stellar center. However, the accurate prediction of 



periods for asteroseismology (for example) must use whole- 
star models so as to account for higher order effects. 

The first several eigenfrequencies are located and 
stored, including the real and imaginary components ujr and 
LOi, the period H and the number of nodes in the eigensolu- 
tion. So far, it has not been necessary to analyze more than 
the first ten eigenfrequencies; the maximum number of un- 
stable modes located so far was nine. The results have been 
collated as diagrams showing the stability of one or more 
modes and the number of unstable modes as a function of 
Toff and L/M. Stability is established by the sign of uJi (neg- 
ative for unstable modes). Results for varying X and Z are 
compared in a mosaic of similar such diagrams, on which 
are superposed the loci of selected groups of variable star. 

When identifying eigenfrequencies it is usual to identify 
the fundamental, first overtone, second overtone, etc, in se- 
quence for each model envelope, as verified by the number of 
nodes in the eigensolution. However, in the most luminous 
models, strange mode pulsations are excited and the one- 
to-one correspondence between the ordinal number of the 
eigensolution identified and the number of its radial nodes 
frequently breaks down. The phy sics of this question are 
discu ssed in more detail elsewhere iSaio. Baker fc Gautschvl 
Il998h : it has no effect on the results presented here. 

Z variation 

The first experiment was carried out on the basis that in- 
creasing the overall metallicity is known to increase the 
size of the Fe-bump instability. Calculations were carried 
out for Z = 0.01,0.02,0.05,0.10 and 0.20. X was reduced 
to 0.80 when X + Z > 1 would otherwise have resulted. 
This is essential l y a re peat of the experiment described by 
ICharpinet et al.l J200j) but over a much larger volume of 
parameter space. 

Figures |21 and |3 show how the resulting instability 
boundaries behave as functions of X and Z for a mass of 
0.50 Mq. The results for other masses (0.4, 0.7 and 0.9 Mq) 
are similar. The discrete model grid is reflected in the con- 
tour shapes. Note that all the panels have hydrogen abun- 
dance increasing from bottom to top and metallicity (in 
whatever form) increasing from left to right. 

For clarity. Fig. |21 shows the instability boundaries for 
only the radial fundamental and the first, second and third 
radial overtones {k = 0,1,2,3). It demonstrates a number 
of important general features. 

1) The Fe-bump instability finger is not significant for 
normal abundances {X = 0.75, Z — 0.02). 

2) The size of the instability finger increases with increasing 
Z and with decreasing X. 

3) For increasing Z, the number of excited nodes increases. 

4) The Fe-bump instability region, including the blue edge, 
becomes bluer for higher order modes (a similar result holds 
for the classical Hell instability strip). 

5) The blue edge for any given mode {e.g. the fundamental) 
is not shifted significantly by increasing Z, while the red 
edge is shifted red- ward. 

6) The low-luminosity end of the instability finger is 
populated by the highest number of excited modes. 

7) Unstable modes with L/M > 3.5 are primarily strange 
modes. 
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Figure 2. Instability boundaries for radial pulsation modes in stars with M = 0.5 Mq and homogeneous envelopes with Z = 0.01 — 0.20. 
The bold line shows the instability boundary for the longest period radial mode, normally the fundamental fc = 0. Thin contours show 
successive boundaries for A: = 1, 2, 3, where these exist. Other features as in Fig. [Tl 



Figure |3] presents an alternative view of the same re- 
sults but, because it includes all the excited modes up to 
k = 10, it reflects the extent and mode-density of the insta- 
bility region more effectively. Many of thes e features were 
already well known from p revious studies jjefferv fc Said 
Il998l: Icharpinet et aljEool . but it is helpful to reiterate 
them here in one unified picture. 

In terms of matching the observations of the EC 14026 
stars and objects like LSIV— 14°116, this figure shows that 



at these high T^s, radial modes (and hence p-modes in 
ge neral) are not excited , even for very high Z (cf. Fig. 3 
in Icharpinet et alj |200l') . Despite the excitation of higher- 
order modes with slightly bluer blue edges, the extension to 
the instability zone is insufflcient. 

This has a natural explanation. In attempting to shift 
the instability boundary, Z has been increased in order to 
increase the driving by the Fe-opacity bump. At the same 
time the mean molecular weight of the gas and the mean 
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Figure 3. The number of unstable radial pulsation modes in stars with M = 0.5 Mq and homogeneous envelopes with X = 0.0 — 0.9 
and Z = 0.01 — 0.20. The bold line shows the overall instability boundary; within this contour at least one radial mode is excited. Each 
subsequent contour (thin lines) represents one additional unstable mode. Other features as in Fig. 



opacity in the layers above the driving zone has been in- 
creased. Therefore the Fe-bump will be located correspond- 
ingly closer to the stellar surface. In order to make the blue 
edge bluer, it is necessary to maintain the contribution to 
driving from iron without increasing the mean molecular 
weight and opacity in the layers above the driving layer. 



Fe variation 

The logical alternative is therefore to increase the iron 
abundance without increasing the abundances of other 
heavy elements. Using the OPAL opacity table generator 
(.Igles ias & Rogers 1996), we constructed new sets of opac- 
ity tables in which the iron abundance alone was increased 
by factors of / = 5, 10 and 20, whilst the remaining ele- 
ments had abundances corresponding to a Z = 0.02 mix- 
ture. There is a modest effect on the overall Z such that 
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Figure 4. As Fig.Elbut with Z = 0.02 and Fe eniianced by factors of 1, 5, 10 and 20. 



Z = 0.02 + (/ - 1)Zfo where Zfc = 0.02188 represents the 
mass fraction of iron within th e standard heavy element mix- 
ture iGrevesse fc Noelslll993h . A similar set of models was 
run for these mixtures. These are compared with the refer- 
ence set Z = 0.02 in Figs.HandU 

It is immediately apparent that increasing the iron 
abundance alone by a factor five produces a predictable ef- 
fect similar to increasing Z from 0.02 to 0.10, but with the 
major difference that the blue edge has become bluer and 
more vertical (in the L / M — T^s diagram). As the iron abun- 
dance is further increased (c/. X — 0.75, Z = 0.02, / = 10), 



a region develops which contains several excited modes and 
which is bluer than the Fe-bump instability finger seen when 
the entire Z-component varies together. 

This is p r ecisely the phenomenon identified by 
ICharpinet et al.l lll997fl . The overall extent of our instability 
zone is larger than identified before; the earher calculations 
benefit from obeying a total iron conservation law by using 
an iron distribution obtained from a self-consistent calcula- 
tion of radiative levitation. 

The significance here is that earlier results are repro- 
duced in character, if not in detail, using a simple model. It 
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Figure 5. As Fig.Elbut with Z = 0.02 and Fe eniianced by factors of 1, 5, 10 and 20. 



tells us that the principal characteristics of p-mode instabil- 
ity in EC 14026 stars may be derived from a homogeneous 
stellar envelope with a factor 10 enhancement in iron abun- 
dance (assuming Z ~ 0.02). This provides useful pointers 
for a subsequent investigation of g-mode instability. 

It is interesting that reducing hydrogen concentration 
extends the instability finger towards lower LjM and lower 
Toff. Whether this region of parameter space is populated by 
any real stars is a matter for conjecture. If they do, they are 
likely helium- and iron-rich white dwarfs with periods of ^ 
seconds. However, diffusion theory suggests that in stars of 



such high gravity, iron would settle out of the layers where 
it would be required to excite pulsations. 

To understand the phenomenon we recall that the in- 
stability blue edge occurs when the thermal timescale 

Tth = CvTAm/L (1) 

becomes too small compared with the pulsation period, 
where T is the temperature in the excitation layer, Cv is 
the specific heat per unit mass at constant volum e, and Am 
is the mass lying above the excitation layer (see ICo:jll980l . 
Ch. 10 for details). Let us assume an ideal gas and that the 
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Figure 7. Fig.lSlredrawn with the observed periods of EC14026 
stars over-plotted on the theoretical periods of unstable radial 
modes for M = 0.50Mq, logL/M = 1.8, X = 0.7b, Z = 0.02 and 
enhanced iron (/ = 10). The modes are labeled by radial degree. 



opacity in the layers above the excitation zone is expressed 
as 



arf-i — 
Kop T 



(2) 



with constants kq, a(0 < a < 1), and b. Then, we can in- 
tegrate approximately, from the stellar surface to the exci- 
tation zone, the hydrostatic equation and the equation of 
diffusive flow of radiation to obtain 



Am oc (KoL^"M")"T+^ri?^ 



(3) 



where is the mean molecular weight and R is the stel- 
lar radius. Substituting equation into equation Q, and 
taking into account Cv oc we have 



Tth CC (KQfJ, 



l + 2a\ 



{L/M)TT^T- 



(4) 



The strong dependence of rth on Tcff is responsible for the 
blue and the red edges of instability range. Since the k- 
mechanism is optimal when the period of pulsation is com- 
parable to Tth, the instability range of a short (long) period 
mode is hotter (cooler) as can be seen in Figs. H and m 

For a given L/M ratio, the blue-edge can be shifted blue- 
ward by increasing rth, i.e. by reducing kq and the mean 
molecular weight /x. The former can be achieved by reducing 
Z without significantly reducing the iron abundance. The 
mean molecular weight is reduced by increasing X. The ef- 
fect is somewhat complicated because an increase in X tends 
to increase kq- Figs. |21 and |1] seem to indicate that the effect 
of reducing the mean molecular weight exceeds the effect of 
increasing opacity. 

Periods 

Predicted pulsation periods for unstable fundamental ra- 
dial modes, and for first, second and third radial harmon- 



ics are shown in Fig. |S|for a range of compositions and for 
logL/M —1,2 and 3. The central column represents mod- 
els with a standard metallicity and reduced hydrogen abun- 
dance, while the left and right columns represent a uniform 
metal enhancement and a selective iron enhancement re- 
spectively. Note again that the lower L/M pulsators for a 
selective iron enhancement are bluer than for a uniform en- 
hancement, and also that higher-order modes are preferen- 
tially excited in bluer pulsators. 

For comparison, the top centre panel shows the pul- 
sation periods and period ranges for several EC14026 and 
extreme helium-star pulsators. A consequence of the models 
is that bluer EC14026 stars should show higher harmonics 
and, conversely, there should be fewer hot sdBVS pulsating 
in the radial fundamental mode. This is clearly refiected by 
the mean slope of the logP(logrcft) relation for EC14026 
stars (excluding PG 1605-1-072 and KPD 1930-f 2752) being 
steeper than that for a single mode. Fig. Qcompares the ob- 
servations for EC14026 variables with pulsation periods for 
unstable radial modes with M = 0.50 M©, X = 0.75, Z = 
0.02 and enhanced iron (/ — 10); logL/M = 1.8 was chosen 
to be representative of the EC14026 variables (see Fig.0. It 
is evident that only the coolest stars are likely to pulsate in 
the radial fundamental mode while many of the hottest can 
only excite the second or third harmonic. Some EC14026 
stars lie beyond the blue-edge for all of these models; a 
higher value of / might have been more appropriate, but 
the overall result is unaffected. 

Two exceptions are PG 1605-1-072 and KPD 1930-1-2752. 
The first has evolved away from the extreme horizontal- 
branch and has a higher L/M ratio (Fig. 0. The sec- 
ond is tidally disto rted by a very close binary companion 
iBilleres et aLllioOCt) and our classical radial models may be 
inadequate. 



Growth Rates 

An indication of the overall instability is given in Fig. |H| 
where the exponential growth times are shown as contour 
plots for compositions representative of extreme helium and 
EC14026 variables. The growth times are represented by 
the dimensionless quantity —cur/i^Oi. Divide by 27r to find 
the growth time in pulsation cycles, and then multiply 
by pulsation period P to obtain the physical growth time 
Tg. Thus, for the EC14026 stars with M = 0.50 Mq and 
logL/M ~ 1.7, Fig. |S1 indicates logP/s ~ 2 and Fig. El in- 
dicates log{—LUr/ijJi) ~ 4, yielding logr^/s ~ 5. For the ex- 
treme helium stars, the same exercise also gives logTg/s ~ 5. 
With growth times ~ one day, these modes are highly un- 
stable. 



3 NON-RADIAL MODES 

Zero-age extreme horizontal-branch stars 

The non-radial modes most usually observed in stars may 
be divided into pressure or p-modes and gravity or g-modes. 
Pressure modes are trapped in the outer stellar envelope and 
may be thought of as a three dimensional generalization of 
the radial modes discussed in the previous section. There- 
fore their stability has essentially been discussed already. 
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Figure 6. Periods of unstable modes due to Fe-bump instability for stars with M = 0.50 Mq and log L/M = 1,2 and 3. Each panel 
shows the fundamental (solid), first (dashed), second (dot-dashed) and third (dotted) harmonics. The middle column shows models with 
a standard metal content {Z = 0.02) and with X increasing upwards (there are no unstable modes for X = 0.75, Z = 0.02 at the L/M 
ratios shown). The left hand column shows models with enhanced Z = 0.20, the right hand column shows models with enhanced iron 
oiily (/ = 10)- Both columns show models with X = 0.0,0.25,0.75 increasing upwards. The loci of selected pulsating stars are shown 
in the centre - top panel; multi-periodic stars are shown by their longest period (square) and their range (line). The selection includes 
pulsating sdB stars and extreme helium stars. 



although the precise periods of p-modes with spherical de- 
gree I > will differ from the corresponding radial (/ = 0) 
values. 

Gravity modes propagate most strongly in the deep in- 
terior of the star. Consequently complete stellar models in- 
cluding both core and envelope are required. Following pre- 



vious work jFontaine et aljr2003l) we have assumed that the 
PG1716 variables are extreme horizontal branch stars and 
have thus constructed series of "zero-age" horizontal-branch 
models having a helium-burning core with a mass of Mc and 
with a hydrogen-rich envelope with mass Me ranging from 
0.0003 to 0.0034 Mq. We have considered two cases of M^; 
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Figure 8. Contour plots representing minimum growth times for unstable modes in stars with M = 0.50 Mq and for compositions 
representative of extreme helium stars (left) and EC14026 variables (right). Contours are separated by 1 dex, the bold lines representing 
log{—ujr/i^i) = 3 and 6, with growth times decreasing toward higher L/M. Other features as in Fig.lTl 



0.476 and 0.486 Mq, The surface layers are characterized by 
mass fractions of hydrogen X and base metalhcity Zq, with 
iron again augmented by a factor /. We have considered se- 
quences with X = 0.90,0.75 and 0.5, with Zq = 0.01 and 
0.02, and with / = 1, 10 and 20. The hydrogen abundance 
Xti at the interface between a hydrogen-rich envelope and 
a helium core is assumed to be given as 

Xtr(Afr) = 0.5X(1 + sine) (5) 

with 

e = -!v[{Mr ~ Mc)/AM (6) 

where Mr is the mass coordinate within the transition layer, 
and AM, which is assumed to be 2 x 10~* M©, is the mass 
of the transition layer. 

The high-temperature end of this sequence corresponds 
with the locus of the short-period EC14026 variables. 
PG1716 variables are to be found at lower temperatures, 
with 25 000 < Tcff/K < 29 000. 

We have carried out a stability analysis for each of these 
models testing the stability of non-radial p- and g-modes 
spherical degree I — 1, . . . , 4. The frequency range consid- 
ered is 0.2 < uj < 20, where u is the angular frequency of 
pulsation normalized by i/ GM/R^ with G being the gravi- 
tational constant. The stability of non-radial modes is found 
to be insensitive to Mc. We show results for Mc = 0.486 M©. 
To allow for the likelihood that most horizontal branch stars 
have evolved away from their zero-age structure, we will 
demonstrate later that their principal stability character- 
istics survive through to the end of core helium burning. 

As anticipated, we find no unstable p-modes with I > 0, 
for any models with normal iron abundance (/ — 1). Like 
the radial modes, these are only excited with significant iron 



enhancements. Unexpectedly, we do find excited g-modes 
/ > 1 in ZAHB stars with Z = 0.02 and no iron enhancement 
over a smaU interval 16 000 < %s /K < 20 000, with I = 2 
modes being excited in higher radial-orders and at lower T^s 
than I = 3 and I = 4 modes (Fig. |5J . 

Some properties of an excited g-mode of Z = 2 are shown 
in Fig. 1101 where the solid line in the bottom panel is the ra- 
dial displacement of pulsation divided by the distance from 
the center ^r/r, the dotted line is the work W, and the top 
panel shows the distribution of the kinetic energy of the 
mode. If a pulsation mode is excited (unstable), W > 
at the surface. There is an excitation zone {dW/dr > 0) 
in a temperature range of 5.1 < log r< 5.4 due to the Fe- 
opacity bump and a damping zone {dW/dr < 0) just below. 
This g-mode is marginally unstable because the driving ef- 
fect is not very strong with no enhancement of iron abun- 
dance. In such a case, a reflection of the wave at the interface 
between the core and hydrogen-rich envelope plays an im- 
portant role. Fig. 1101 shows that the kinetic energy and the 
amplitude of the mode suddenly decrease towards the cen- 
ter at log T ~ 7.25, where the hydrogen abundance changes 
steeply, indicating a partial reflection of the mode to occur 
there. A sudden change in the mean molecular weight at the 
interface requires a sudden change in the gas density, which 
reflects g-modes partially. The reflection reduces the ampli- 
tude in the core and hence reduces radiative damping there. 
The strength of the reflection depends on the phase of the 
spatial oscillation at the interface, and hence only some of 
the g-modes are excited in the optimal period range for the 
Fe-bump mechanism. 

One model for which g-modes are excited on the ZAHB 
was evolved through its complete horizontal-branch evo- 
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Figure 9. Periods of modes due to Fe-bump instability for ZAHB stars with A/c = 0.486 M©, X = 0.75 and Zo = 0.02, and Mrme = 
increasing from right to left. The left hand column shows models with no additional iron enhancement (/ = 1), the right hand column 
shows models with iron increased by a factor / = 10. Stable modes are marked as (blue) dots, unstable modes are marked as filled (red) 
circles. The chemical mixture with / = 1 corresponds to that in the first column, second row of Figs. |^ and |^ while that with / = 10 
corresponds to that in the third column, second row of Figs. Blandly 



lution. The g-modes continue to be excited throughout 
HB evolution, primarily because, with no hydrogen-burning 
shell, the basic structure of the envelope and the character- 
istics of the core-shell reflection surface are not altered. 

The above result indicates that g-mode pulsations can 
be driven in certain blue horizontal branch stars irrespec- 
tive of whether chemical stratification enhances the iron- 



abundance in the driving zone. This conjecture will require 
observational verification, which will constrain theoretical 
models for the blue horizontal branch stars. Moreover, it im- 
plies that g-modes in these cooler stars are more likely to be 
excited with modest iron enhancements than the p-modes. 
The growth times for the ZAHB models are relatively long 
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Figure 11. As Fig. El for models witii Mc = 0.486 Mq, X = 0.75 and iron enliancement / = 20, but for Z = 0.02 (left) and O.Ol(right), 
demonstrating the shift of the g-mode blue edge with metallicity Z. 



(~ 100 y, Fig. llUL but are reduced by factors of ~ 10 during 
HB evolution. 

However, in general, unstable modes (both p- and g-) 
are only recovered with a significant enhancement of iron. 
Indeed, with / ~ 10, we always find unstable modes of ei- 
ther p- or g-type down to T^s ~ 13 000K (Fig. El . In these 
models, the refiection at the core boundary is not important 
in the excitation of g-modes, and most modes are excited 
in an optimal range of periods. The growth times for these 
modes are ^ 1 — 10 y. The locus of the p-mode instability 



follows that of the radial modes as a function of X, Z and 

In our models, with / > 10, the high temperature 
limit of the unstable g-modes always overlaps the low- 
temperature limit of the unstable p-modes of the same 
spherical degree. While this overlap is marginal for I — 1 
modes, it increases substantially with I. This result has an 
important corollary which states that if a star has an ef- 
fective temperature such that it is in the p- and g-mode 
instability overlap (p/g overlap) and if p-mode oscillations 



Fe-bump instability: the excitation of pulsations in subdwarf B and other low-mass stars 13 



X = 0.90, Z = 0.01, f = 20 



X=0.50, Z = 0.01, f = 20 



o 

• I — I 

0) 

o 



3 4 

o 

. I — I 

O 



O 

. 1=1 

Oh 

o 




J I I I I I I I 

4.3 4.2 
log 



Figure 12. As Fig. El for models witii Mc = 0.486 Mq, Z = 0.02 and iron enliancement / = 20, but for X = 0.90 (left) and 0.50(right), 
demonstrating the shift of the g-mode blue edge with hydrogen fraction X. 



are observed, i.e. it is an EC14026 variable, then the iron 
enhancement necessary to drive the p-modes will also be 
sufficient to drive the g-modes, i.e. it is also likely to be a 
PG1716 variable. 

This is true regardless of the fact that we have used a 
global iron enhancement rather than a stratified model. The 
reason is that the driving for both p- and g-modes comes 
from the same part of the star, so that for a star of the 
right effective temperature, it may actually be impossible 
to have p-mode instability without g-mode instability. In- 



deed, the stars Balloon 090100001 and HS 0702-1-6043 lie on 
the boundary between the observed EC14026 and PG1716 
variables and exhibit both types of pulsational behaviour 
JSchuh et alJl2006l: lOreiro et alJl2005ll . 

However, there is a problem. The high temperature 
limit of the unstable g-modes occurs around Teff ~ 20 000 
K for I = 1 modes, increasing to ~ 25 000 K for I = 4. 
The observed high temperature limit for PG1716 variables 
is ~ 29 000 K and it is unlikely that I > 2. 

We have investigated g-model stability for ZAHB mod- 
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Figure 10. The bottom panel shows the radial displacement t;,. 
and work W of an unstable high-order g-mode of ( = 2, and 
opacity k as a function of the temperature. In a driving (damping) 
zone, dW/dr > (< 0), and the mode is unstable if > at 
the surface. The top panel shows the distribution of the kinetic 
energy of the mode. A sudden change in the kinetic energy and in 
the amplitude of the spatial oscillation at logT 7.25, where the 
hydrogen abundance changes steeply. The convective core, where 
the g-modes are (spatially) evanescent, extends to log T = 7.934. 
G.R. is the growth rate —uii/LUr- 



els having a range of chemical composition. Figs. Ill| and ll2l 
demonstrate the effects of varying the overall metallicity Z 
or the overall hydrogen concentration X whilst maintaining 
a high iron concentration (/ = 10). Increasing X (equiv- 
alent to reducing the helium fraction) has the effect of a 
modest increase in the effective temperature of the g-mode 
blue edge. Similarly, reducing Z while conserving the iron 
abundance (Fig. IllH shifts the blue edge of g-modes blue- 
ward a little. These effects can be understood from equation 
101 However, in neither case is a substantial change in the 
envelope composition capable of shifting the g-mode blue- 
edge to its observed location. 

Figs. initoIiS exclude resuffs for Z = 4. The I = 4 blue 
edge is always bluer than the I = 3 blue edge by ~ 0.01 
dex. In addition to the lower likelihood of observing I = 4 
modes, this is again insufficient to account for the observed 
blue edge. 

These results clo sely mirror those obtained previously 
l)Fontaine et al.ll2003fl . where radiative levitation provides 
the necessary enhancement in the iron abundance in the 
driving region. These authors found g-mode instability with 

1 = 3,4 at effective temperatures up to 24 000 K (log Tea < 
4.38), or with Z < 8 up to logTeff < 4.43, stiU too cool 
to account for the observed blue edge. The main difference 
is that our models show unstable g-modes ai I — 1 and 

2 where those of Fontaine et al. did not. One reason may 
be that our model envelopes have a homogeneous chemical 



Figure 13. Evolution tracks for horizontal branch stars with 
Mc = 0.476 M0, X = 0.75, Zq = 0.02 and enhanced iron abun- 
dance / = 10. Total masses are M = 0.4763, 0.478, 0.480, 0.483 
and 0.490 Mq. Models having either unstable p- or g-modes with 
I < 4 are indicated by (a) and (▼), respectively. Small dots indi- 
cate models in which neither is excited. Note that models having 
both unstable p- and g-modes appear as "stars". 



composition with an artificial iron enhancement, whereas 
the earlier models have a stratified composition with iron 
enhanced in critical layers by a well-understood mechanism. 

We have found that the effective temperature at the 
blue edge of g-modes hardly changes when the degree of 
iron enrichment increases as long as / > 10. This indicates 
that to have a blue edge of g-modes consistent with the 
observation, we have to change something in the structure 
of models. 

Evolved horizontal-branch stars 

One solution to the blue-edge question might be to con- 
sider the evolution of these stars. We have computed the 
evolution of iron-enriched extended-horizontal branch mod- 
els with X = 0.75, Zo = 0.02, / = 10 and Mc = 0.476 Mq 
through the completion of core helium-burning and beyond 
(Fig. 1131 . Note the "breathing pulses" (small loops in the 
HR diagram) towards the end of core helium-burning. These 
are encountered in our models with artificially high iron 
abundances. They are due to the increased core opacity and 
are not significant in models with a normal core composi- 
tion. S e micon vection is included following the description by 
ISpruid lll992l) 

We have carried out a stability analysis for selected 
models along each evolutionary sequence. Fig. 1131 shows 
those models where either p-modes and/or g-modes are 
unstable for I < 4. Comparing the instability tempera- 
ture range of zero-age models for the same composition. 



Fe-bump instability: the excitation of pulsations in subdwarf B and other low-mass stars 15 



we see that stability of p-modes are hardly affected by 
the evolutionary change of the stellar structure, while g- 
modes tend to be damped in evolved models in which the 
central part is radiative. There is a clear overlap region 
4.30 < logTeff/K < 4.38, where both p- and g-modes are un- 
stable. This extends upwards in luminosity to include some 
post-horizontal-branch stars; it does not extend blue-wards. 

The Brunt- Vaisala frequency becomes extremely large 
in the evolved non-degenerate radiative core. This makes the 
wavelength of a g-mode very short, which in turn enhances 
the radiative damping and stabilizes the mode. Still some 
g-modes are excited even in evolved models with a radiative 
core. For these models a partial reflection at the bottom of 
the helium-rich envelope reduces the amplitude and hence 
radiative damping in the core. Since only few g-modes are 
excited and the evolutionary speed is fast, these g-modes are 
not very important observationally. 



4 CONCLUSION 

Our original goals were to understand oscillations observed 
in the helium-rich subdwarf LS IV— 14° 116, and in the 
PG1716 variables. We therefore carried out a broad review 
of Fe-bump driven pulsational instability in low-mass stars^ . 

By first considering radial modes, we have demon- 
strated the essential role of chemical composition such that 
instability increases with the contrast between the iron- 
bump opacity and other opacity sources. This increased con- 
trast may be achieved either by increasing the iron ab un- 
dance, confirming earlier work bv lCharp inet et aL] l|200l|) or 
by reducing the hydrogen abundance dJeffe^^^Sai d ll99d) . 
At least one of these is necessary to excite oscillations in all 
of the low-mass Fe-bump pulsators discovered to date. 

We have further demonstrated that the blue-edge for 
radial instability is affected by the mean molecular weight 
in the stellar envelope, so that increasing the iron abundance 
alone provides a bluer instability region than increasing the 
iron abundance in concert with all elements heavier than 
helium. The former is required to explain the locus of the 
EC14026 instability region, making the general assumption 
that the properties of non-radial p-modes modes are closely 
linked to the corresponding radial mode of the same radial 
order. 

Furthermore, the blue-edge also depends on the radial 
order of the oscillations, so that higher-order modes may be 
found in hotter stars. By comparing theoretical and observed 
periods for EC14026 stars, we have shown that low-order or 
fundamental mode radial oscillations are only likely to be 
seen in the coolest EC14026 stars and that the hottest stars 
must oscillate in higher-order modes. 

However, we have been unable to explain the oscillations 
in the helium-rich star LSIV— 14°116; the observed periods 
are simply too long compared with the L/M ratio obtained 
from spectroscopy. A lower surface gravity would deflnitely 
help. 

Considering non-radial pulsations, we have focused on 
g-mode instability and in particular on instability in modes 

^ Oscillations in main-sequence stars were deliberately excluded 
from our study. 



with spherical degree I < 4. We have discovered a small g- 
mode instability island on the blue horizontal branch which 
does not require iron enhancement or hydrogen depletion. 
However, since it requires Z — 0.02 and a partial reflection of 
the wave at the interface between the hydrogen-rich envelope 
and the helium core, it will take some observational effort 
to verify whether there are any real pulsating horizontal- 
branch stars which correspond with these models. Effective 
temperatures are expected to be between 16 000 and 20 000 
K, and pulsation periods between 1 and 4 hours. 

With a factor of 10 enhancement of iron on a back- 
ground metallicity Z = 0.02, a large instability region de- 
velops, even for 1 = 1, which extends from < 13 000 to 
~ 24 000 K (Z = 3). With such iron enhancement there is 
always an overlap between p- and g-mode instability regions 
so that stars near the boundary would be expected to exhibit 
both modes simultaneously. Depleting the envelope hydro- 
gen abundance tends to shift both the g-mode blue edge and 
the radial/p-mode red edge to lower temperatures. 

Many of these results have been demonstrated individ- 
ually ('Saio'1993'; 'Je fferv fc SaiJl998l : ICharpinet et al .119961 : 
fFontainc et al. 200I|); this investigation has placed them in a 
more general context, as well as delivering some new results. 
Of these, the role that mean molecular weight plays in de- 
termining the instability boundaries for radial and p-mode 
oscillations has implications for understanding the oscilla- 
tions in EC14026 variables. We have made no effort to jus- 
tify the adopted abundances on physical grounds, we have 
simply sought parametric solutio ns which satisfy the obser- 
vatio ns. Radiative levitation ( Clia ver. Fontaine fc WesemaeJ 
[1995") is known to operate in extreme horizontal-branch stars 
and provides a natural explanation for the required iron en- 
hancements. The fact that it operates selectively accords 
well with our deduction that only iron should be enhanced. 

Our discovery that g-mode oscillations may be excited 
in blue horizontal branch stars with envelopes of "normal" 
composition has two-fold consequences. In addition to the 
observational question already posed, it may be less hard 
to excite g-mode pulsations in PG1716 stars than hitherto 
supposed; maybe we don't have to work so hard to find the 
necessary chemical structure in the stellar envelope as we do 
in the case of EC14026 variables. 

Nonetheless, we have been unable to resolve the prob- 
lem of why the observed boundary between the EC 14026 
and PG1716 stars occurs at ~ 29 000 K. The theoretical p- 
mode/g-mode boundary remains persistently at ~ 24000 K 
over a wide range of envelope compositions. It is possible 
that a sharp discontinuity in composition immediately be- 
low the iron bump could help; this and other experiments 
with stratified envelopes remain to be investigated. 
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